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Abstract 

The problem of wave propagation in a biaxial graded-index fiber with circular symmetry is 
considered. The problem is formulated in terms of four first-order differential equations for 
the tangential components of the electric and magnetic fields. A general solution method 
for solving systems of differential equations is presented. This solution method is then used 
to solve the system of equations for a particular example of a biaxial graded-index fiber. 
Numerical results for the propagation constant in the fiber are also given. 



I. Introduction 

The optical fiber has become a much studied transmission system due to its property of 
wave guidance with low loss. In recent years it has been shown that introducing anisotropies 
into the dielectric medium of the fiber produces several interesting features, such as control 
of power flow and reduction of peak attenuation near cutoflf. 

Typically the analysis of wave propagation in a cylindrical dielectric waveguide such 
as an optical fiber is performed using a wave equation formulation. For the simple case 
of a step-index fiber a detailed analysis, including dispersion relations, cutoff conditions 
and mode designations, is presented by Snitzer [1]. Paul and Shevgaonkar [2] present a 

similar analysis for a uniaxial step-index fiber and also perform a perturbation analysis to 
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deter min e the modal attenuation constants. These are the only two cases for which exact 
solutions are known. 

For inhomogeneous fibers no exact solutions are known. For the case of an isotropic 
graded-index fiber several approximate analytic solution methods are available. These ap- 
proximate solutions all share the common assumption that the fiber is infinite in extent. 
In addition if the permittivity is assumed to vary slowly over the distance of one wave- 
length the wave equation formulation simplifies to an associated scalar wave equation. If 
the permittivity profile is parabolic the solution to the scalar wave equation can be written 
in terms of either Laguerre polynomials [3] if cylindrical coordinates are used or Hermite 
polynomials [4] if rectangular coordinates are used. For arbitrary permittivity profiles the 
scalar wave equation can be solved using the well known WKB solution method [5], [6]. For 
a parabolic permittivity profile all three solution methods give identical results. Under the 
assumption that the fields are far from cutoff Kurtz and Streifer [7], [8] have shown that a 
solution to the full vector problem can be written in terms of either Laguerre polynomials if 
the permittivity profile is quadratic or asymptotically in terms of Bessel and Airy functions 
for arbitrary permittivity profiles which decrease slowly and monotonically. A comparison 
of the vector and scalar solutions for the quadratic permittivity profile implies the vector 
modes can be obtained by simply renumbering the scalar modes [9]. Usmg the renumbered 
scalar modes as a basis Hashimoto [10] and Ikuno [11] have developed two slightly differ- 
ent iterative methods which can be used to solve the full vector problem for an isotropic 
graded-index fiber. 

An alternate formulation of the problem is to write the four first-order differential 
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equations for the tangential field components as a first-order matrix differential equation. 
For a step-index fiber with uniaxial core and cladding Tonning [12] has shown that the 
matrix formulation can be solved exactly in terms of Bessel functions. For isotropic graded- 
index fibers with arbitrary permittivity profiles Yeh and Lingren [13] have indirectly used 
the matrix formulation in developing a numerical solution method based on the concept 
of stratification. Using the concept of transition matrices Tonning [14] has developed a 
numerical procedure which can be used to solve the matrix differential equation for isotropic 
graded-index fibers. 

II. Formulation of the Problem 

Consider a circularly symmetric optical fiber with the geometry shown in Figure 1. In the 
core, 0 < p < a, the permittivity is given by e 0 € r {p) where to is the permittivity of free space 
and e T {p) is the relative permittivity tensor of the core and is a function of p only. In the 
cladding, a < p < b, the permittivity is given by e 0 e c where e c is the relative permittivity 
of the cladding and is assumed to be constant. In both the core and the cladding the 
permeability is p 0 , the permeability of free space. For convenience, the external radius of 
the cladding, b , is assumed to be sufficiently large in comparison to the radius of the core, 
a, so that it is not necessary to impose boundary conditions at the air-cladding boundary. 
Consider the case where the relative permittivity tensor in the core is given by 

e r (p) = 0 e 2 {p) 0 ; ( 1 ) 

\ 0 0 e 3 (p)/ p ,<M 

where d{p), e 2 (/>) and t 3 {p) are the relative permittivities in the p, <f> and 2 directions 
respectively. In general the relative permittivities are arbitrary functions of p. However, 
the choice of cylindrical coordinates requires ei(p) and e 2 (p) be equal at p = 0. 
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For time harmonic fields in a source free region, Maxwell’s equations can be written as 

V x H = ju>e o^rE, ( 2a ) 

V x E = -jvp oH, ( 2 &) 

where w is the angular frequency. If the z and 4> dependence of the fields is given by 

where p is the longitudinal wavenumber and m is any integer, then for cylindrical coordi- 
nates Maxwell’s equations can be written in component form as 


— H z + = u>€ot\Epi 

9 

(3a) 

dH z 

-jPH p - - JW€ 0 €2^, 

(36) 

l±( p E4,)- 3 —H p = j»e 0 e 3 E z , 

pdp v p 

(3c) 

— E z + fiE# = -<jjp 0 H p , 

P 

(3d) 

dE z 

jPE p + = jupo H+, 

(3e) 

l±(pE<,)- 3 —E p =-ju>PoH z . 

( 3 /) 


The remainder of the problem can now be formulated in two different ways. If the transverse 
field components E p , £<*, H p and H* are eliminated from eqs. (3) we obtain a pair of 
coupled second-order differential equations for the longitudinal field components E z and 
H z . Alternately, if the radial components E p and H p are eliminated we obtain a system of 
four first-order differential equations for the tangential field components E z , E+, H z and 

H#. 
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First consider the coupled wave equation formulation. It is convenient to define a 
normalized magnetic field h = Z 0 E where Z„ = %/W^ « the impedance of free space. 
Solving eqs. (3a), (3b), (3d) and (3e) for E p , Ej,, h p and gives 



(4a) 

(46) 

(4c) 

(id) 


where k 0 = is the free space wavenumber and k\ n = k%<i n {p) - /3 2 , n = 1,2 is 

the transverse wave number. Substituting the expressions for E p , E$, h p and given by 
eqs. (4) into eqs. (2c) and (2f) and making a change of variable from p to a normalized 
radius r = p/ a results in the following pair of coupled differential equations for E z and h t 


E" + f\(r)E z + A 2 gi(r)E t = p 2 (r)h' z + q 2 (r)h z , (5a) 

K + M r )K + A 2 Ja(r)fr, = Pi(r)E' z + qi(r)E z , (56) 


where ' — d/dr, A 2 = (fco a ) 2 i K — P/ko ^d 


.. x 1 * V l( r ) 

/l(r) 'r C i(r)[ €l (r)-^], 

(6a) 

... 1 

h( r ) — r ( 2 (r) — k 2 

(66) 

cfrl. . . „L m 2 € 2 (r) ] 

ffl(r) ” €,(r) €l r * ] [ A 2 e 3 (r)[e 2 (r) - **]r» J, 

(6c) 

9i(r) = h(r) -«][!" A 2[ Cl ( P ) _ K 2y\ 

( M) 

, v jmK [fi(r) - e 2 (r) 

pi(n- r [ €l(r) _ K 2j 

(6c) 
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P i(0 
?l( r ) 
92(1 r ) 


jmK 

[<? 2 (r) - €i(r) 

(l{r)r 

r 

On 

to 

1 

* 

to 

jmK 

r 

«1W 1 

T 

[e 2 (r) - k 2 J ( 

jmn 


«iM 1 

e a (r)r 

* 

1 

fH 

w 
1 


( 6 /) 

(69) 

( 6 /») 


The equations for E z and h z become uncoupled for three particular cases. For the so-called 


meridional modes m is equal to zero and therefore from eqs. (6e/,g,h) so are the functions 
Pl (r), p 2 ( r ), qi (r) and g 2 (r). For isotropic and uniaxial step-index fibers e x and e 2 are equal 
and constant and again from eqs. (6e,f,g,h) the functions pi(r), p 2 {r), ?i(r) and q 2 {r) are 


zero. 


In general a solution of eqs. (5) for arbitrary permittivity profiles is not possible. It 
is possible to obtain a fourth-order differential equation for either E z or h z by eliminating 
h z or E z from eqs. (5). However, the complexity of the resulting equation precludes the 
determination of a solution. For meridional modes a direct series solution of the uncoupled 
equations is possible. However, due to the poles in the functions / x (r) and / 2 (r) the resulting 
series solution will not be convergent for the entire core region. An exact solution of eqs. (5) 
is possible only for the case of a step-index fiber. For either an isotropic or uniaxial step- 
index fiber the coupled equations simplify to Bessel’s differential equation. 

In order to find an analytic solution of eqs. (5) some assumptions must be made. First, 
the cladding is neglected and the core is assumed to extend to infinity. This eliminates the 
need to impose boundary conditions on the solution at the core-cladding boundary. Second, 
the permittivities are assumed to be slowly varying functions of r over a distance of several 
wavelengths. This is equivalent to assuming e-(r) ft 0. For the case of either an isotropic or 
a uniaxial graded-index fiber, e x (r) = c 2 (r), application of the second assumption to eqs. (5) 
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results in the following equations for E z and h z 

E';+-E' z + A 2 gi{r)E z = 0, ( 7a ) 

T 

K + ~K + A 2 g 2 (r)h z = 0, (™) 

T 



(8a) 

( 86 ) 


For the case of a biaxial graded-index fiber, e a (r) / c 2 (r), the previous assumption does 
not cause eqs. (5) to uncouple since pi(r) and p 2 (r) not identically equal to zero. 

Eqs. (7) can be solved easily using the well known Wentzel-Kramers-Brillouin (WKB) 
solution method [5], [6], [15]. The solutions obtained using the WKB method are not 
solutions of the full vector problem given by eqs. (5) but rather they are solutions of a 
related scalar problem given by eqs. (7). However, the vector solutions can be obtained by 

renumbering the solutions to the scalar problem [9]. 

For the case of a biaxial graded-index fiber the WKB solution method can be applied 
blindly to eqs. (5) and at most only two terms in the WKB expansion can be determined. 
Nevertheless, the term representing the phase of the WKB solution is not a well behaved 
function and therefore it is not reasonable to assume the WKB solution method remains 
valid under this condition. In order to solve Maxwell’s equations for the case of a biaxial 

graded-index fiber an alternate formulation must be used. 

Instead of eliminating the transverse field components from eqs. (3), eliminate the radial 
field components E p and H p and write the remaining four equations as a system of four 
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first-order differential equations in terms of the tangential components [14]. Prom the two 
algebraic equations, eqs. (3a) and (3b), the radial components can be written in terms of 
the tangential components as 


E p = 
B P = 


1 

u>e 0 €i 

1 

U>/io 


jE M +fiH 4 ] 


(9a) 

(96) 


Using eqs. (9), the four remaining equations can be written as 

iSl = -j™h, + -L(„ - «’)(,)*). (10a) 

ds S€ 1 5€i 

= — (™ 2 - fi« 2 )^ + j^(*M, ( 106 ) 

^irzj^^-^ez-K 2 )^) ( 10C ) 

ds s S 

(10J) 

(15 5 5 

where a change of variable from p to a normalized radius s = k 0 p has been made. Eqs. (10) 
can be written in matrix form as 

$5 = - A(s)u, (H«) 

ds s 


where 


>2 

u = ( E z sE^ h z shj , ) 


(1U) 


and 


/ 


A (*) = 


0 

0 


-JTT 


J-(m 2 - «i s 2 ) 


jmK -j{( 2 ~ « 2 ) 

\-j{m 2 - ( 3 s 2 ) -jmK 


0 

0 


0 

0 


(11c) 


For the special case of meridional modes, m = 0, eqs. (10) can be separated into two 


systems each containing two equations. The first set corresponding to transverse magnetic 
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modes (TM) can be written in matrix form as 


du (TM) = 1 A (TM) (3)u (TM) 


ds S 


(12a) 


where 


u(™> = (E t sh, 


(126) 


and 


A<™>(*) = ( ° 

\ je 3« 


0 £(«i - « 2 ) 


2 0 


(12c) 


The second set corresponding to transverse electric modes (TE) can be written as 

du — = Ia( te >(s)u< te > 

ds s 


(13a) 


where 


u< TE > = (/i, «£*)’ 


(136) 


and 

«=(->.= ' i( T !) ). (13C) 

The only known exact solutions of the matrix equation are for the cases of an isotropic and 
a uniaxial step-index fiber [12], [14]. These solutions are identical to the exact solutions of 
the wave equation formulation. 

It is not readily apparent that the matrix equation is easier to solve than the wave 
equation formulation. As was mentioned earlier, a series solution for the wave equation 
formulation is possible only when the equations are uncoupled. However, for the meridional 
modes of a graded-index fiber no series solution will be convergent for the entire core region. 
In contrast, the system matrix A(s) does not have any poles in the core region and therefore 
the series solutions will be convergent in the entire core region. 
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While the form of A(s) guarantees a convergent series solution the series may not con- 
verge rapidly enough to use it in numerical computations. An alternate solution method is 
asymptotic partitioning of systems of equations [16]. This method involves the transforma- 
tion of a system of linear first-order differential equations into a system of equations whose 
solutions are easier to find. The form of the solution method presented in the next section 
is based on the expansion of the a general system matrix A(x) in terms of positive powers 
of x, in contrast to the usual form where the expansion is in terms of powers of l/x (see 

e.g. [16]) 

III. Matrix Partitioning 

Consider the following system of N linear differential equations 

— = — A(a:)u(i) as x — > 0 (14) 

dx r 9 

where u is a column vector, q is an integer greater than or equal to 1 and A(x) is a N x N 
matrix given by 

A(i) = 52 AnXU “ x °- 

n=0 

It is possible to simplify this system of equations by transforming them into some special 
differential equations whose solutions are easier to find. Let 

u(z) = P(*)v(*) ( 16 ) 

where v is a col umn vector and P(x) is a N xN nonsingular matrix. Using eq. (16), the 
original problem given by eq. (14) can be transformed into 


( 17 ) 



where 
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B(x) = P(x)' 1 A(x)P(x) - (18) 

or more conveniently 

= A(x)P(x) - P(x)B(x). ( 19 ) 

ax 

The matrix P(x) is chosen so that B(x) has a convenient form, either diagonal or Jordan 
canonical form. If B(«) has either of these forms the solution of the transformed system 
for v is trivially obtained. For example, if A(x) is a constant matrix then P(z) is also a 
constant matrix and eq. (18) is simply a similarity transformation. This implies P(*) is 
chosen so that B(x) is either the diagonal or the Jordan canonical form of A(x). In general, 
when A(x) is not a constant matrix P(x) is not a constant matrix and it is not dear from 
either eq. (18) or (19) how P(x) should be chosen so B(x) has the desired form. 

■ In order to develop a procedure to find B(x) and P(x) start by expanding them as the 
following Taylor series 

OO 

B(x) = B n x” as x -> 0, 

n =o ( 20 ) 

OO 

P(*) = £ p „*" as x — > 0, 

n=0 

where in general B 0 is a Jordan canonical matrix and B„ is a diagonal matrix. Substitutmg 
eqs. (20) into eq. (19) and equating like powers of x gives 

A 0 Po - PoBo = 0 ( 21 ) 

for x° and 

n 

(ri — q + l)P n _g+i = ^ (A/P n «/ - PjB n _/) 

l=o 


( 22 ) 
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for *», n > 1, where P n ., +1 = 0 for n - q + 1 < 0. Eqs. (21) and (22) define an iterative 
procedure to find the coefficient matrices for the series expansions of B(i) and P(z) so that 
either eq. (18) or (19) is satisfied. Eq. (21) can be rewritten as 

Bo = P^AoPo ( 23 ) 

which implies P 0 is chosen so that B 0 is either the diagonal or Jordan canonical form of 
A 0 . With some algebraic manipulations eq. (22) can be written more conveniently as 

B 0 W n - W n B 0 = (n - q + l)W n _,+i + B n - F n (24) 

where the matrices W n and F n are defined as 

W n = Po ( 25 ) 

and 

F n = PJ 1 A„P0 + P.-' E (An— (P| ' P.B„_,) . (26) 

/ = l 

Notice that the unknowns in eq. (24) are the matrices B n and W n and that the matrices 
W n _, + i and F n depend solely on matrices found in previous iterations. Since by definition 
B n is a diagonal matrix eq. (24) can be solved easily for B„ and W n by setting the diagonal 
elements of B n and F n equal to each other and then solving for W n from what remains of 
eq. (24). Since the form of B 0 is known in advance an explicit solution in terms of B 0 and 
F n can be found for W n 

Consider the special case where B 0 is a diagonal matrix and q = 1. This corresponds 
to the form of the matrix differential equation (11) which we want to solve. While it is not 
obvious from eq. (11c) that B 0 will be a diagonal matrix, it will be shown later that this is 


indeed the case. 
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When B 0 is a diagonal matrix the expression B 0 W„ - W„B 0 has zeroes along it s 
main diagonal and does not depend upon the elements along the main diagonal of W n . 
The solution to eq. (24) can be easily written as 


(B1 i = '’ ; 

( n)i> ‘ \ 0, i # J, 


(27) 


and 


i = j; 


(28) 


(Wn)ij i * j > 
where A;, t = 1, 2, . . . , N are the eigenvalues of A 0 . One potential problem exists with this 
solution. If Ai - A j - n = 0 and (F n ) t , # 0 for some particular values of i, j and n then it 
may not be possible to find W n and therefore a solution may not be possible. 

Consider a biaxial graded-index fiber with permittivity profiles of the power law type 


given by 


Cj(r) = €j(l — 2A;r“’) i = 1,2,3 


(29) 


where ei = 6,(0) and Ai - (f, - 6 c )/2e,. Since the choice of the coordinate system reqmres 
fi(0) = e 2 (0), the definition of A; requires that Aj and A 2 be equal. Then, for this choice 
of permittivity profiles e,(r) and e 2 (r) are not equal only when oo # a 2 . The case of a 
step-index fiber exists as a special case to the power law profiles in the limit as a, - oo, or 
equivalently by setting A; = 0. 

Now let us solve the matrix equation for the transverse modes in a biaxial fiber where 
the permittivity profiles are parabolic. The relative permittivity profiles can be written in 
terms of the normalized radius s = kop = (k 0 a)r as 


tj(s) = 6{(1 - 2A°s 2 ) i — 1,2,3 


(30) 



14 


where A? = A t /(fc 0 a) 2 . Strictly speaking this choice for e;(s) does not produce a biaxial 
fiber since by definition cj = e 2 and A? = A®. Since e x (s) and «,(*) do not appear together 
in the matrix equations for the transverse modes, eqs. (12) and (13), it is not necessary to 
set e 2 equal to e a . However, in the final result it is necessary to replace c 2 by e x and either 
set A x = 0 and obtain the solution for a biaxial fiber where e x («) is constant and « 2 (s) is 
parabolic or set A° = 0 and obtain the solution for the case where c x (a) is parabolic and 
f2 (s) is constant. The first term in the series expansions for A (TM) («) and A (TE) (<) is given 


by 



(31a) 


wVl PTf* 

i = ™ ; (316) 

1 -j(€ 2 -K 2 ), 1 = TE. 

Since the two eigenvalues of A 0 are both equal to zero it is not possible to find a non-zero 
P 0 such that Bo is a diagonal matrix. Instead P 0 must be chosen so that B 0 is a Jordan 
canonical matrix. For A c as given by eq. (31) choose P 0 as 


so that 



(32) 


(33) 


is a Jordan canonical matrix. After four iterations, the solution for the TM case is found 


to be: 
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+ 


sh,*, = —Ci < 


a 2 [Vci> 


c iVl 


/lz = ~j*N 2 C l 

„ r k 2 

sE * = -Ci 


> 64 €1 3 8 J J 

(34a) 

16 €i 3 2 J J 

(346) 

The solution for the TE case is 


*N 2 .2 , k N2 g 2 
4 64 


(35a) 

_ *£2^^ 4 


(356) 


2 16 

where = 62 - k 2 . 

Now consider the solution of the matrix equation for hybrid modes. For all permittivity 
profiles the first term in the series expansion of A(s) is 


An = 


/ 0 0 

0 0 

jmK ~jk 2 m 
\-jm 2 -jmK 


k 2 v 

- m* j Z£LL \ 

J <i 
-^2 

€l 

0 

0 0 / 


0 


•TFIIt 

Jir 


(36) 


where Jk 2 ^ = e a - k 2 and the eigenvalues of A 0 are ±m, m # 0. Since the eigenvalues are 
repeated, in general, the choice for P 0 should at best cause B 0 to be a Jordan canonical 
matrix. This is the only restriction placed on the form of P 0 by the solution method. Any 
P 0 which causes B 0 to be a Jordan canonical matrix can be expected to result in a valid 
solution. Since it is possible for several different choices of Po to satisfy this condition, 
conceivably there may exist several possible mathematical solutions to the problem. 

Since the solution for a step-index fiber exists as a special case of the solution for a 
graded-index fiber it is reasonable to choose P 0 based on the knowledge of the exact solution 
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for a step-index fiber. Prom the wave equation formulation we know that for a step-index 
fiber the differential equations for E z and h x become uncoupled and the resulting equations 
can be solved independently of each other. This suggests that for the case of a step-index 
fiber P(s) and hence P 0 should have a form such that two of the four elements in the 
solution of the vector v(s) should contribute to E t but not h x while the remaining two 
elements contribute to h z only. If Po is chosen as 


Po = 


( ^Nl 

0 


0 \ 

mK 

jm 

mK 

-jm 

0 

^Ni 

0 


V -jme 1 

mK 

jme i 

mK / 


(37) 


then for a step-index fiber E z and h z are at least uncoupled for the lowest order solution 

where P(s) = Po- 
ll sing P 0 given by eq (37) B 0 is given by 


Bo 


(m 00 0 \ 

0 m 0 0 

0 0 — m 0 

V 0 0 0 -m ) 


(38) 


Since Bo is a diagonal matrix, instead of a Jordan canonical matrix as was the case for the 
transverse modes, eq. (28) can be used to find W n . Recall that this solution for W„ may 
cause some elements of W n to be undefined. In particular, for this problem the elements 
in the third and fourth columns of both W 2 jt and P2fc undefined when m — 1,2, ... ,k. 
However, due to the structure of the various matrices and the order of multiplication in the 
definitions of W„ and F„ these undefined elements remain in the third and fourth columns 
of all resulting matrices. In the final solution these undefined elements can be dropped since 
they contribute only to the two solutions which are not finite at s = 0. 

With P 0 and B 0 given by eqs. (37) and (38) respectively, the general form of the solution 
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to eq. (11a) which is finite at s = 0 is given by 


E t \ 


{ Pu{ a ) Pu{ s )\ 

sEj, 

— c m 

P 2 l{ a ) P 2 i{ a ) 

h z 

— 5 

P3\{*) Ptti*) 

\ sh 4 ) 


\Pa i{») Pa2(s)J 


V C 2 e A »(*> 


) 


where 


(39a) 


N 

Paw = £(*.)«*" 

n=0 

" s n 

M 5 ) = 53 ^ n ^**"n 

n— 1 

and N is the number of iterations. 

For a biaxial graded-index fiber where ei(s) and e 3 (s) have a parabolic profile given by 
eq. (30) and e 2 (s) is a constant, after two iterations the following expressions are found for 

A i(s) and Pij(s) 


i = 1,2, 3, 4; J = 1, 2 

(396) 

i = 1,2 

(39c) 



(40a) 


(406) 

= (<1 _ «-) + 4rii(t ^ + 1} [;> («, - «=)= - ^ . 

(40c) 


(40(f) 

**<*>- ”"' + 4(m + l){«i ( * I_K ^ 

+ "^* A ?h- «*) + (« + «*.]}**• 

(40e) 


+ [("* + l)* 2 - ( f i - * 2 )] } a2 > W) 
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Psii*) = ~j 


■ m /o aO\ .2 


(2A?)5 2 , 


4 (m + 1) 

= (<. - * ! ) + 4m( ^ +1) [('> - « 2 ) ! - • 
P 41 (,) = -jmei + - « ! ) - - *"_ + K ‘ ),e 

P„(„) = m* + [(«i - *’) - ”* |” + j )fl (2A;)]» J • 


(40j) 
(40*) 
(40i) 
(40j ) 


This should not be considered an accurate solution for u(s) since the term A° does not 
appear anywhere in eqs. (40). This solution is identical to the solution obtained after two 
iterations for a biaxial graded-index fiber where €i(a) has a parabolic profile and e 2 (s) and 
f 3 (s) are constant. Since A§ only appears in the matrix A 4 at least four iterations must be 
performed in order to obtain the effects of a non-constant f3(«)- 

The solution for a uniaxial or a step index-fiber can be obtained from eqs. (40) by 
setting A?( and A°) equal to zero. Notice that setting A? equal to zero causes P J2 (s) 
and P 31 (s) to be set equal to zero. This corresponds to the decoupling of the differential 
equations which occurs in the wave equation formulation for the case of a step-index fiber. 

From n um erical results, it appears that the functions A i(s) and Pij(s) given in eqs. (40) 
are monotonic. This indicates that the solutions for the various field components will not 
have an oscillatory behavior. Consequently, for a given value of m only the mode with the 
lowest cutoff frequency will be found. 
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IV. Numerical Results 

As was previously stated the solution for the biaxial graded-index fiber given by eqs. (40) 
does not include the effects of a non-constant e 3 («). Obtaining a more accurate solution 
requires performin g more than four iterations. Instead of deriving algebraic equations for the 
elements of F„, B n , W n and P„ the values of these matrices can be determined numerically 
if the values of m, k and k 0 a are known in advance. One potential difficulty with this method 
comes from the undefined elements in W n and P n . Since these elements contribute only to 
the solutions which are unbounded at s = 0 they can be set equal to zero without affecting 
the final solution. The ability to do this appears to depend upon the form of A(s) and the 
ordering of the eigenvalues of Ao in Bo. 

Asymptotic partitioning was used to solve the matrix equation for several types of 
fibers. For the case of a step-index fiber a comparison was made between the propagation 
constants determined using asymptotic partitioning and those determined using the exact 
solution. For transverse modes the asymptotic solutions were in poor agreement with the 
exact solutions. Since for transverse modes in a step-index fiber asymptotic partitioning 
produces a series solution, the poor agreement can be attributed to using too few terms 
in the series expansion of the exact solution. For hybrid modes there was a much better 
agreement between the asymptotic solutions and the exact solutions. In particular, for a 
step-index fiber the asymptotic and the exact solutions produced almost identical values 
for the propagation constants of the HEn mode. 

Figure 2, 3 and 4 are plots of the normalized propagation constant for the HE n , HE 2 i 
and HE31 modes in a biaxial graded-index fiber where ei(s) and e 3 (s) have a parabolic 
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profile and f2( a ) is a constant. 

V. Conclusions 

For both the wave equation formulation and the matrix equation, exact solutions are known 
only for the cases of an isotropic step-index fiber and a uniaxial step-index fiber. For 
isotropic and uniaxial graded- index fibers the wave equation formulation can be solved 
approximately using WKB analysis. For a biaxial graded-index fiber WKB analysis can not 
be used on the wave equation formulation. Asymptotic partitioning can be used to solve 
the matrix equation for all types of permittivity profiles. For meridional modes asymptotic 
partinning appears to generate the series solution for the matrix differential equation. For 
hybrid modes, of the solutions produced by asymptotic partitioning have a form such that 
for a given value of m only the mode with the lowest cutoff frequency can be found. A nice 
feature of the asymptotic solutions is that they remain valid all the way down to cutoff. 
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Figure 2 Normalized propagation 
fiber where £i(s) and «3( 
constant . 


constant for the HE n mode in a biaxial graded-index 
s) have the parabolic profile of eq. (30) and «2(-») is a 
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Figure 3 


Normalized propagation constant for the HE 21 mode in a biaxial graded-index 
fiber where €l (a) and t 3 {s) have the parabobc profile of eq. (30) and e 2 (s) is a 

constant . 
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Figure Captions 


Figure 1 Geometry of the fiber 

Figure 2 Normalized propagation constant for the HEn mode in a biaxial graded-index 
fiber where d (a) and € 3 ( 5 ) have the parabolic profile of eq. (30) and t 2 («) is a 
constant. 

Figure 3 Nor mali zed propagation constant for the HE 2 i mode in a biaxial graded-index 
fiber where ei(s) and e 3 (s) have the parabolic profile of eq. (30) and e 2 (s) is a 
constant. 

Figure 4 Normalized propagation constant for the HE 31 mode in a biaxial graded-index 
fiber where e a (a) and e 3 (s) have the parabolic profile of eq. (30) and € 2 (s) is a 


constant. 



